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Abstract 

In this article, we study the radiative decays among the charmonium states with the heavy 
I ' quark effective theory, and make predictions for the ratios among the radiative decay widths 

^ of an special multiplet to another multiplet. The predictions can be confronted with the 

. experimental data in the future and put additional constraints in identifying the X, Y, Z 

' charmonium-like mesons. 

^ ; PACS numbers: 14.40.Pq; 13.40.Hq 

. Key Words: Charmonium states, Radiative decays 

<^ : 

CN . 1 Introduction 

(— I , In recent years, the Babar, Belle, CLEO, DO, CDF and FOCUS collaborations have discovered (or 

i' confirmed) a large number of charmonium-like states and revitalized the interest in the spectroscopy 

of the charmonium states. Many possible assignments for those states have been suggested, such as 
^ • the conventional charmonium states, the multiquark states (irrespective of the molecule type and 

<~i , the diquark-antidiquark type), the hybrid states, the baryonium states, the threshold effects, etc 

[mills, 4, 5, GjZi- The main difficulties in identifying those new mesons as the excited charmonium 
, states in that the observed masses do not fit in with the predictions of the potential models based 

^ ' on the confining potential. In this article, we will focus on the traditional charmonium scenario 

, and study the radiative decays among the charmonium states with the heavy quark effective theory 

■ [iisiiin]- 

In 2003, the Belle collaboration observed the X(3872) in the J/^tt+tt" invariant mass distribu- 
tion with a significance in excess of 10 cr in the B if J/-07r+7r~ decay [11], and later the X(3872) 
was confirmed by the DO, CDF and Babar collaborations [HllISlIT^. The decay X(3872) J/^-/ 
' observed by the Belle and Babar collaborations [HI (El [17] and the decay X(3872) — >■ ip'j ob- 

, served by the Babar collaboration |17j favor the even charge conjunction assignment. The di-pion 

r* ' spectrum also indicates that the charge conjunction is even [11] , and the studies of the Belle collab- 

. !^ I oration [TS] and CDF collaboration [T51 [TH] favor the spin-parity-charge-conjunction J^'^ — 1++ 

' assignment for the X(3872). However, the recent analysis of the B — > J/ipujK decay by the 

^ I Babar collaboration indicates that the P-wave orbital angular momentum for the J/ipu system 

■ - - ' is more favored than the 5*- wave, the X{3872) maybe have the spin-parity-charge-conjunction 

J^'-' — 2 ^ instead of the J^*-^ — 1++ [20^. In Ref.|21). Jia et al study the radiative decays 
X(3872) J/^p{t(j') + 7 within several phenomenological potential models, and observe that the 
2 ^ assignment for the X{3872) is highly problematic. In this article, we assume that the X(3872) 
is the conventional charmonium Xci(2P) [22, 23, 24l f25l [26l . 

In 2005, the Belle collaboration observed the ^(3930) in the DD invariant mass distribution 
near 3.93 GeV in the 77 collision with the statistical significance of 5.3 cr [27], the mass and width 
are Mz {3930) = (3929 ± 5 ± 2)MeV and rz(393o) = (29 ± 10 ± 2) MeV, respectively. The pro- 
duction rate and the angular distribution favor the Xc2(2P) assignment |27| . In the same year, 
the Belle collaboration observed the X(3940) in the recoiling spectrum of the J/tp in the process 
e+e^ — > J/V' + D*D [55]. Later the Belle collaboration studied it with higher statistics and de- 
termined the mass and width Mx(394o) = (3942^^ ± 6) MeV and rx(394o) = (37t^| ± 8) MeV 
|29) . Furthermore, they observed the X(4160) in the D*D* invariant mass distribution in the 
process e~^e~ — >■ J/ip + D*D* with a significance of 5.1 a. The mass and width of the X(4160) are 
Mx(4ieQ) = (4156^20 ± 15) MeV and rx(4i60) = {^^^-IV ± 2l)_MeV, respectively. The observa- 
tion of the dominant decay mode of the X(3940) being the D*D and the lack of evidence for the 
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DD decay mode indicate that it is a good candidate for the ?7c(3S) state [35J [511 130] : although the 
mass is lower than the prediction of the potential models [31] . The dominant decay mode of the 
X(4160) is the D*D* and the DD and DD* modes have not been observed, so the X(4160) 
may be the 77c (4S) or Xco(3P) candidate [32] . 

Also in 2005, the Babar collaboration observed the y(4260) in the J/i/^Tr+Tr" invariant mass 
distribution in the initial state radiation (ISR) process e'^e~ — >■ 'yis'RJ/'4'T^~^'^~ [33] . It was con- 
firmed by the CLEO collaboration ^34, and Belle collaboration ^5], and in the Belle data there 
is also a relative broad structure y(4008) with the mass 

My(aoo8) = (4008 ± 40^28^) MeV and 
width ry(4oo8) = (226 ± 44 ± 87) MeV, respectively. In Ref.|36], Llanes-Estrada suggests that the 
y(4260) may be the ^^(48) charmonium state and displaces the -0(4415). 

In 2006, the Babar collaboration observed the structure y(4320) in the initial state radiation 
process e+e~ — 7isr'0'''' with the mass My(432o) = (4324 ± 24) MeV and width ry(432o) = 
(172 ± 33) MeV, respectively [37], and later the Belle collaboration observed two relative nar- 
row resonant structures F(4360) and y(4660) with the masses My(43go) = (4361 ± 9 ± 9) MeV, 
My(4660) = (46 64 ± 11 ±5) MeV and the widths ry(436o) = (74 ± 15 ± 10) MeV, ry(466o) = 
(48 ± 15 ±3) MeV, respectively [3S]. In 2008, the Belle collaboration observed the r(4630) as 
a threshold enhancement in the A+Aj invariant mass distribution in the initial state radia- 
tion process e+er ^ 7isrA+A;: [^. The mass and width are My(463o) = {4.GMtrtl) MeV 
and ry(463o) — (92^24^21) MeV, respectively, which are roughly in agreement with that of the 
y(4660). In Ref.[5^, Li and Chao calculate the mass spectrum of the charmonium states based 
on the screened potential, and observe that the Z(3930) agrees with the Xc2(2P), the ■0(4415) 
is compatible with the V(5S) rather than the i/'(4S), the r(4260), r(4360) and y(4660) may be 
the 0(4S), -0(30) and 0(6S) respectively, and the X(3940) and A(4160) may be the ??c(3S) and 
Xco(3P) respectively. In Ref.[40^, Ding et al identify the y(4360) and r(4660) as the 3^01 and 
5'^Si charmonium states respectively (Badalian et al share the same opinion [41]). and evaluate 
the e~^e~ leptonic widths. El transitions. Ml transitions and the open flavor strong decays. 

In 2009, the CDF collaboration observed a narrow structure F(4140) near the J/00 threshold 
with a statistical significance in excess of 3.8 a in the exclusive B — >■ J /^(j)K decays produced in pp 
collisions [42]. The mass and width are MY{iUQ) = (4143.0 ± 2.9 ± 1.2) MeV and ry(4i4o) = 
(ll.7t^:g±3.7) MeV, respectively [42]. The r(4140) is very similar to the charmonium-like 
state F(3940), which was observed by both the Belle and Babar collaborations near the J/'0w 
threshold in the exclusive B J/ipuK decays [43l [44]. The mass and width are Mr (3940) = 
(3943 ± 11 ± 13) MeV and ry(394o) = (87 ± 22 ± 26) MeV respectively from the Belle collabora- 
tion [ig and Mr(394o) = (3914.6l|j ± 2.0) MeV and ry(394o) = (341^^ ± 5) MeV respectively 
from the Babar collaboration [^. In 2009, the Belle collaboration reported the observation of 
a significant enhancement with the mass (3915 ± 3 ± 2) MeV and total width (17 ± 10 ± 3) MeV 
respectively in the process 77 — > ujJ/tp [15], these values are consistent with that of the y(3940). 
The updated values of the mass (3919.ll|^ ± 2.0) MeV and total width (31^8° ± 5) MeV from 
the Babar collaboration are also consistent with the old ones [46]. The Belle collaboration mea- 
sured the process 77 — > 4>J/ip for the J li^4> invariant mass distributions between the threshold 
and 5 GeV, and observed a narrow peak A"(4350) with a significance of 3.2 cr [47]. The mass and 
width are Mx(435o) = (4350.6l^ j ± 0.7) MeV and rx(435o) = (l3.3l^^i^ ± 4.l) MeV, respectively; 
and no signal for the F(4140) J/ipip structure was observed [47]. It is difficult to identify the 
y(3940), F(4140) and AT (4350) as the conventional charmonium states [48] . 

In Table 1, we list the experimental values of the charmonium states with the possible identifica- 
tions compared with the theoretical predictions [26l [SU [49] . We do not mean that such assignments 
are correct and exclude other possibilities, and just take it for granted for the moment and study 
the radiative decays among the charmonium states with the heavy quark effective theory [8l [H [10] , 
which have been applied to identify the excited Dg and D mesons, such as the Z?s(3040), Ds{2700), 
£1,(2860), L»(2550), £1(2600), £1(2750) and L»(2760) [50l [H] [52l [53] [54] . 

The article is arranged as follows: we study the radiative decays among the charmonium states 
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Table 1: Experimental and theoretical mass spectrum of the charmonium states. The SP, NRP and 
GI denote the screened potential model, the non-relativistic potential model and the Godfrey-Isgur 
relativized potential model, respectively. 
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with the heavy quark effective theory in Sect. 2; in Sect. 3, we present the numerical resuhs and 
discussions; and Sect. 4 is reserved for our conclusions. 



2 The radiative decays with the heavy quark effective the- 
ory 

The charmonium states can be classified according to the notation n^*+^Lj, where the n is the 
radial quantum number, the L is the orbital angular momentum, the s is the spin, and the j 
is the total angular momentum. They have the parity and charge conjugation P = (— 1)-^+^ 
and C = (—1)-^+^, respectively. The states have the same radial quantum number n and orbital 
momentum L can be expressed by the superfields J, J^, J^"^, etc [55] . 



J = 
J'^ = 



l + f 
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1 



V2 



Xa + r]c2 75 



(1) 



where the denotes the four- velocity associated to the superfields, the fields xfn^ai X%j xji^ 
x", have the total angular momentum j — 3,2,1,0,2 respectively, and belong to different 
multiplets. The fields in a definite superfield have the same n, and form a multiplet. We multiply 
the charmonium fields xLa: xL, Xn, X°, Vc2 ^ ' • • with a factor y/M^ 



and they have 



dimension of mass |. The superfields J/^i - z^j- have the following properties under the parity, charge 
conjunction, heavy quark spin transformations. 



p 

C 

s 
p 



7°Jmi...m.7°, 



(-1)^+^C[J,,...^J^C, 
e J e't 



(2) 



where S, S' S SU (2) heavy quark spin symmetry groups, and [S, f] = [S", f] — 0. 

The radiative transitions between the m and n charmonium states can be described by the 
following Lagrangians, 



^ss 
C-sp 

CpD 



E 
E 
E 

m.Ti 



^(m,n)Tr[J(m)(7^,J(n)] 

8{m, n)Tr [J(m) J^(n) + J^(n) J(m)] V^" , 

5(m, n)Tr [J^,(m) J"(n) + J^H J^,(m)] ^ 



(3) 



where J^j...^^ = 7°'^/ti m-l'^^' ~ F^^v^j, the F^'' is the electromagnetic tensor, and the (5(m, n) 
are the coupling constants, which have different values in the Lag rangians C^ss^ ^sp^ ^pd^ 
use the same notation for simplicity. The Lagrangians Csp and Cpo preserve parity, charge 



4 



conjugation, gauge invariance and heavy quark spin symmetry, while the Lagrangian Css violates 
the heavy quark symmetry. There are two tensors (J^vF^''{— a^vFap) and ct^i/F'"' can be 

used to construct the spin-breaking Lagrangian Css ■ They both have the terms a ■ E and a-B, and 
can turn into each other with the interchange E -(r^ B, we choose the parity-conserving structure 
a^vF^^ in constructing the spin-breaking Lagrangian. We take the Lagrangian Csp from Ref.[H] 
and construct the Lagrangians Css and CpD, and later Dr. F. De Fazio draws my attention to 
Ref. [56j , where the Lagrangian Cp^ is introduced for the first time. 

From the heavy quark effective Lagrangians Css, ^SP and Cpo, we can obtain the radiative 
decay widths F, 

F = ^^y-^|T|2, (4) 

where the T denotes the scattering amplitude, the kj is the momentum of the final states in the 
center of mass coordinate, the J2 denotes the sum of all the polarization vectors, the j is the total 
angular momentum of the initial state, and the Mi is the mass of the initial state. Cho and Wise 
study the radiative decays of the heavy quarkonia using the multipole expansion and heavy quark 
symmetry [57j . 



3 Numerical Results 

We calculate the radiative decay widths F using the FeynCalc to carry out the sum of all the polar- 
ization vectors. In calculations, the masses of the charmonium states are taken as the experimental 
values from the Review of Particle Physics [H], see Table 1; for the unobserved charmonium states, 
we take the values from the screened potential model |26[ , and assume that the charmonium states 
V'CS^Si) and 77c(5^So) have degenerate masses. In this article, we identify the i/)(4415) as the V'(5'^Si) 
charmonium state, and expect that the 77c(5"'^So) charmonium state has slightly smaller mass. The 
predictions of the screened potential model are M^(4g) = 4463 MeV and A/,j^(4g) = 4446 MeV, 
respectively, so we take the approximation M^(4S) = M^^(4s) = 4421 MeV. 

The numerical values of the radiative decay widths are presented in Tables 2-4, where we retain 
the unknown coupling constants S(rn,n) among the multiplets of the radial quantum numbers m 
and n. In general, we expect fitting the parameters S{m, n) to the precise experimental data, 
however, in the present time the experimental data are far from enough. In Tables 4-6, we present 
the ratios of the radiative decay widths among the charmonium states. 

The radiative decay widths listed in the Review of Particle Physics are 
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(5) 
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where we have neglected the uncertainties [491. From those radiative decay widths, we can obtain 
the following ratios, 



r(V>(2S) -> Xco(lP)7) 
r(7^(2S)-^Xci(lP)7) 
r(7/;(2S)-^Xco(lP)7) 
r(^(2S) ^ Xc2(lP)7) 

r(xco(iP) ^ J/^7) 
r(xc2(iP) ^ J/^7) 
r(xci(iP)^ J/^7) 
r(xc2(iP) ^ J/^7) 

r(V^(3770) Xco(lP)7) 
r(V'(3770)^Xci(lP)7) 
r(V>(3770)^Xco(lP)7) 
r(V'(3770)^Xc2(lP)7) 



^ 1.046(1.151), 

= 1.101 (1.649), 

= 0.311 (0.363), 

= 0.770(0.756), 

= 2.517(3.167), 

> 8.111 (82.02), (6) 



where the values in the bracket are the theoretical calculations based on the heavy quark effective 
theory. The agreements between the experimental data and the theoretical calculations are rather 
good, and the heavy quark effective theory works rather well. The ratios presented in Tables 4-6 
can be confronted with the experimental data in the future at the BESIII, KEK-B, RHIC, PANDA 
and LHCb, and put additional constraints in identifying the X, Y, Z charmonium-like mesons. 
There is a relative P-wave between the final-state charmonium and the photon, the radiative 

decay widths T oc k'^, where fc-y = — ^m- * ' ^^'^ ^'^f denote the masses of the initial and final 

charmonium states respectively. The numerical values of the decay widths shown in Tables 2-4, 
where the effective coupling constants 6(m,n) have been factorized out, reflect the corresponding 
processes are facilitated or suppressed in the phase-space. If the energy gap Mi — Mf is small 
(or large), small variations of the masses Mi and Mf can (or cannot) lead to remarkable changes 
for the decay width. For example, we plot the radiative decay widths r(Ar(3872) J/4'l) a-nd 
r(A:(3872) V'7) versus the mass Mx in Fig.l. 

From the experimental data of the Babar collaboration Br(S+ A:(3872)if +) xBr(X(3872) 
J/iPj) = (2.8 ± 0.8 ± 0.2) X 10-^ and Br{B+ X{3872)K+) x Br(A:(3872) ^ V'7) = (9.9 ± 2.9 ± 
0.6) X 10~^ [17], we can obtain the central value of the ratio ^rfx^3872)^i^-y) ~ 0-283, which means 
^2(2 2) ^ 0.00576. The large hierarchy 5(2, 1) <^ S{2, 2) is compatible with the phenomenological 
expectation that the Xci(2P) state potentially decays to the tp'j rather than to the J/ipj, 

The El and Mi transitions among the charmonium states are usually calculated by the formula 

[SilMllMlini], 



|2 



^ / 2s+l r . /2s' + l^, \ 4 2 a 2j' + 1 ,2s' + 1, I 2s+l r \ |2 

Tmi (n + ^ n ^ / + 7 ) - '^M^2LTT '"'"^^ ' ^fl^ 



(7) 



where the Ej is the photon energy, the Ef is the energy of final state charmonium, the Mi is the 
mass of the initial state charmonium, and the angular matrix factor Cfi is 

Q, = max(L,L')(2/ + l)| ^' J }' • (8) 

The values of the matrix elements (n'^" ^^L'jr \ r \ n^^^^Lj) and (n'^" ^^L'ji \ n^^'^^Lj) depend on 
the details of the wave-functions which are evaluated using a special potential model, for example. 
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Table 2: The radiative decay widths of the S'-wave and D-wave charmonium states to the P-wave 
charmonium states. The unit is 10^^S^{m,n). 

the Cornell potential model, the logarithmic potential model, the power-law potential model, the 
QCD- motivated potential model [62], the relativized Godfrey- Isgur model, the non-relativistic 
potential model [3T], the screened potential model [13], the relativistic quark model based on 
a quasipotcntial approach in QCD [63l [64], etc. All predictions should be confronted with the 
experimental data. In this article, we intend to make estimations based on the heavy quark 
effective theory, and prefer the ratios among the radiative decay widths of an special multiplet to 
another multiplet, where the unknown parameters J(m, n) are canceled out with each other. 

4 Conclusion 

In this article, we study the radiative decays among the charmonium states with the heavy quark 
effective theory, and make predictions for ratios among the radiative decay widths of an special 
multiplet to another multiplet, where the unknown couple constants 5(m, n) are canceled out with 
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Table 3: The radiative decay widths of the P-wave charmonium states to the S'-wave charmonium 
states. The unit is 10~'^ S^{m,n). 
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Table 4: The ratios of the radiative decay widths of the S'-wave charmonium states to the S'-wave 
charmonium states. The unit of the widths is 5^(m, n). 





r(V' ^ Xci7) 


r(V' ^ Xco7) 




2S^ IP 


1.433 


1.649 


1.154 


2S-)- IP 


1.053 


1.101 




3S^ IP 


0.766 


0.390 


1.195 


3S^2P 


2.042 


1.087 


0.056 


4S^ IP 


0.708 


0.311 


1.895 


4S 2P 


0.934 


0.379 


1.918 


4S-^3P 


2.176 


1.428 


3.090 


5S^ IP 


0.688 


0.283 


1.931 


5S^2P 


0.810 


0.308 


2.006 


5S->3P 


0.876 


0.373 


2.443 


6S^ IP 


0.669 


0.256 


1.607 


6S ^ 2P 


0.732 


0.263 


1.490 


6S^3P 


0.716 


0.267 


1.343 


ID -s- IP 


25.91 


82.02 




ID -s- IP 


> 3.222 


> 8.111 




2D ^ IP 


18.77 


35.08 




2D ^ 2P 


28.85 


50.81 




3D ^ IP 


18.15 


30.54 




3D 2P 


21.40 


33.17 




3D -> 3P 


25.27 


46.71 





Table 5: The ratios of the radiative decay widths of the 5- wave and D-wave charmonium states 
to the P-wave charmonium states. There we normaUze F^^-^^^-y = 1) the wide-hat denotes the 
experiential values. 
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r(xci ^ il^i) 


r(xco ^i) 




IP^ IS 


0.756 


0.363 


1.549 


IP^ IS 


0.770 


0.311 




2P^ IS 


0.842 


0.744 


1.262 


2P -> 2S 


0.462 


0.279 


1.360 


3P^ IS 


0.952 


0.866 


1.148 


3P ^ 2S 


0.854 


0.748 


1.124 


3P ^ 3S 


0.566 


0.342 


2.803 



Table 6: The ratios of the radiative decay widths of the P-wave charmonium states to the S- 
wave charmonium states. There we normaUze T-^^^^^-^ = 1, the wide-hat denotes the experiential 
values. 




Figure 1: The radiative decay widths of the X(3872) versus its mass, the A and B denote the 
processes X(3872) — J/V'T and X(3872) — > i/''7, respectively; the units of the decay widths are 
10-2^2(2^1) and 10-3^2(2,2), respectively, (5(2,2) > (5(2,1). 
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each other. The predictions can be confronted with the experimental data in the future at the 
BESIII, KEK-B, RHIC, PANDA and LHCb, and put additional constraints in identifying the X, 
Y , Z charmonium-like mesons. 
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